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1 Introduction 



Multicritical matrix models have served as an important tool in the study of 2d 
quantum gravity coupled to matter. They were introduced by Kazakov for this 
purpose [1] and the identification of the first multicritical matrix model limit as 
corresponding to a (p, q) = (2, 5) minimal rational conformal field theory coupled 
to 2d Euclidean quantum gravity was due to Staudacher [2]. 

The concept of multicritical so-called branched polymers (BP) was introduced 
in an attempt to understand the physics underlying the multicritical matrix mod- 
els [3]. Ordinary branched polymers were encountered in the theory of random 
surfaces which aimed at providing a non-perturbative definition of the Polyakov 
path integral, either on a hyper-cubic lattice [5] or using the formalism of dynam- 
ical triangulations (DT) [6\. For these non-perturbatively defined random surface 
theories the result was that in physical dimensions (D > 2) the random surfaces 
degenerate into branched polymers [5j [7] . The same phenomenon was observed 
in attempts to study higher dimensional Euclidean quantum gravity using DT 
[5] . Branched polymers seem to be very generic structures which are entropically 
favored and in statistically inspired models of random geometry one has to make 
special efforts to avoid them. 

The formalism of causal dynamical triangulations (CDT) represents an at- 
tempt to eradicate the dominance of BP. It has been partly successful as a model 
of higher dimensional quantum gravity [9]. The two-dimensional CDT model [TU] 
is special in the sense that it is exactly solvable and in the sense that it can be 
understood as a specific limit of ordinary Euclidean quantum gravity where baby 
universes have been integrated out [11]. It is somewhat amusing that while one 
of the purposes of the CDT model was to avoid BP, the two-dimensional model is 
entirely described by BP, as first noticed in [T2]. We can use the BP description 
of two-dimensional CDT in a constructive way to formulate multicritical CDT 
models. In the following we will do that, and we will show that these multicritical 
CDT models are special cases of more general multicritical matrix models. The 
scaling limit taken for these matrix models, the so-called classical scaling limit, is 
different from the conventional scaling limit referred to above. This new scaling 
limit was first introduced in a matrix model with the purpose of obtaining CDT 
from a matrix model [15] . We present here the generalization to multicritical 
behavior. 
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Figure 1: The graphical representation of eqs. ([7]). The black colored graphs 
represent the class of BPs where the link touching the root contains no dimer. 
The red-colored graphs represent the class of BPs where a dimer touches the root. 
It is assumed that only one link is incident on the root vertex. This restriction 
can easily be lifted, but is imposed as it leads to slightly simpler algebra. 

2 Multicritical branched polymers and multi- 
critical CDT 



2.1 Multicritical branched polymers 

One can define a statistical ensemble of BP by the partition function 



(i) 



BP i 



The sum is over graphs which are rooted connected planar trees. The root is a 
distinguished vertex. For simplicity we assume there is only one link incident on 
the root. This assumption has no consequences for the critical behavior of the 
ensemble of BP. To each vertex i is assigned a weight Vi which we usually take to 
depend only on the order of the vertex. To each link / is assigned a fugacity e _M . 

The graphical selfconsistent equation which determines Z{p) is shown as the 
lower part of Fig. [TJ with the colors ignored, in the special case where v\ = v% = 1, 
all other being zero. 

It leads to the following equation 



1 + v 2 Z + v 3 Z 2 



f(Z) 



F(Z). 



(2) 
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We assume by simple rescaling that v\ = 1. The generic BP is obtained if 
the weights Vi are non-negative^] In this case F(Z) has an minimum Z where 
F'(Zq) = and F"(Zq) > 0. Thus we have the non-analytic behavior 

Z(//)-Z(// )~(/i-/io) 1/2 , e^ = F(Z ), (3) 

in the neighborhood of the critical point /io- 

However, if we give up the requirement that the weights Vi of branching should 
be positive it is clear that there exists special choices of the Vi such that not only 
F'(Z ) = 0, but also F^(Z ) = 0, k = 1, . . . m - 1. In this case we obtain 

Z<^)-Z<jj*)~c<^-iM>) x l m , e^ = F(Z ). (4) 

For m > 2 we say that the ensemble of BP is multicritical [31 [EH 0]. The 
fractal structure of the BP depends on m. The Hausdorff dimension of a type m 
multicritical BP is 

m 

d h (m) = -. (5) 

m — 1 

As m increases the BP gets more and more one-dimensional at the critical point. 

The multicritical behavior of the BP can be given a concrete realization as 
hard dimer models with negative fugacities on BPs with positive weight, very 
much like the situation for the multicritical matrix models, as pointed point in 
|15j . We define a BP hard dimer model as 

^(M) = £11^)11^ E £' HD(BP)I - ( 6 ) 

BP i I HD(BP) 

Here we will assume all non-negative and for each BP we sum over all hard 
dimer configurations. The dimers live on links in the BP, "hard" meaning that 
the dimers are not allow to touch. To each dimer we associate a fugacity £. For a 
given dimer configuration HD the total weight will thus be ^ HD ^, \HD\ denoting 
the number of dimers in the configuration HD. To illustrate this let us consider 
the simplest BP model with positive weight, namely i>i,t>3 = 1, all other weights 
equal zero. Thus F(Z) = (1 + Z 2 )/Z and the model indeed has a critical point, 
Hq = log(2), Z = 1. Now allow hard dimers on the links, and let the fugacity of 
a dimer be £. As shown in Fig. [T] we have now two rooted partition functions, 
one starting with a link without a dimer £)), the other starting with a link 
with a dimer (Z(/i, £)). The corresponding equations are: 



1 One can obtain non-generic behavior for non- negative weights Vi if infinite many > 
and F(Z), Z > 0, is a monotonic decreasing function in the range where it is defined 17 . 
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For positive fixed values of £ one has the standard critical behavior of Z as a 
function of /i, i.e. a critical point /io(£) such that near that point (|3| is valid. 
However, writing pi) as 



V - MO = 02(e) - ZMO)) 2 + <*(£) - zOn,®)) 8 + ■ ■ ■ , (8) 

the coefficient C2 will decrease as £ becomes negative, and eventually one will 
reach a point £ c , with a critical point /i c = //o(Cc) where C2(£ c ) = 0. At the 
corresponding Z c = Z (/i c ,£ c ) we have 



<9Z 



(9) 



dZ 2 

Clearly this point is a m = 3 multicritical point: 

li- /i c ~ c 3 (Z(n,£ c ) - Z c ) 3 , i.e. Z((jL,£ c ) - Z c ~ c(/i - /i c ) 1/3 . (10) 

One can calculate the critical exponent related to dimers, i.e. the matter 
system. Following [HE] the critical //o(0 has the interpretation as the free energy 
density of the matter system and the critical exponent o (corresponding to the 
magnetization at the critical temperature in the case of a magnetic system) can 
be defined as 

^ (11) 

where the term on the rhs denotes the leading non-analytic term when £ ap- 
proaches £ c . It is straight forward to show that in addition to (10) we have 

£-£ c ~(Z(// (0,0-^) 2 , (12) 

from which we deduce that a = 1/2. 

The multicritical behavior reported above in terms of dimers is generic for BPs 
and independent of the particular model. Another explicit BP model, relevant for 
the theory of causal dynamical triangulations, is one where branching of arbitrary 
high order, with weight one, is allowed. The graphical equations are shown in 
Fig. [2] and the equations corresponding to (|7]) are: 

e "4(rV(T^)- ^=«rb- (13 > 

Also this model can be solved explicitly and the critical behavior is as above. 

The simplest dimer generalization to higher multicritical points is obtained by 
considering several kinds of dimers with independent fugacities and allow them 
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Figure 2: The graphical representation of the equations (13). 



to touch the same vertex, but not to be put on the same link. One can then 
arrange for the fugacities such that one obtains a m = 4 multicritical BP. 

We note the following: an Ising model coupled to a generic BP cannot be 
critical |16j . Although the Hausdorff dimension of the BP is two, part of the 
inherent structure of the BP is still too one-dimensional to allow for a critical 
behavior. Thus Ising spins cannot change the fractal properties of an ensemble 
of BP. However the above calculation shows that coupling hard dimers to BPs 
can lead to a change of the critical properties of branched polymers. We will 
comment on the criticical properties of the dimers below. 



2.2 Multicritical CDT 



The ordinary two-dimensional CDT theory is defined by the following partition 
function 



'CDT 



W = E' 



(14) 



T 



where the summation is over the set of so-called causal dynamical triangulations 
T, usually defined with two boundaries, separated a link distance t. We might 
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Figure 3: The construction of a tree on a CDT configuration. The inner circle is 
the entrance boundary. The interior to the inner circle is just an additional con- 
struction to obtain bijection between the ensemble of BPs and the CDT ensemble 
of random surfaces (see [12] for a discussion). 

contract the boundaries to contain only a point if convenient. For detail we 
refer to [1 Oj . There is a bijection between causal triangulations and BPs, as first 
observed in [12]. Let us mention how this assignment is done. To each vertex in a 
CDT triangulation, except those at the two boundaries we have a number of links 
pointing forward in time, a number of links pointing backwards in time and two 
spatial links. The assignment is now that all links pointing forward from a vertex, 
except one link which one can take as the link to the left, belong to the BP. For 
the entrance boundary, which has one marked vertex, a special assignment has 
to be made. For details we refer to [12] and for an illustration see Fig. [3j 

This BP reaches all the vertices of CDT triangulation and to each CDT trian- 
gulation there is a distinct BP. Similarly to each BP one can construct a distinct 
CDT triangulation. The partition functions of the ensembles of BP and CDT are 
thus identical except for a trivial redefinition of /i. 

From this identification it is clear that one can add a matter system on 
the CDT triangulation which changes the critical behavior of the triangulation, 
namely one can add hard dimers on the BP associated with the CDT triangula- 
tion. Adding such dimer systems to the CDT ensemble will change the Hausdorff 
dimension of the CDT ensemble at the multicrititical dimer point, such that a 
m-multicritical dimer system, m = 3,4, . . . leads to a CDT random surface sys- 
tem with Hausdorff dimension (pj). That the Hausdorff dimension of the CDT 
ensemble is larger than or equal to that of the BP ensemble is clear, since there 
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are more link-paths connecting two points in a CDT triangulation than in the 
corresponding the BP-ensemble. If we assume there is a unique fractal structure 
of the CDT ensemble at all scales at the critical point, it also has to agree with the 
fractal dimension assigned to the BP. This is because the shortest link distance 
from the root of a BP to any vertex is essentially the same as in the corresponding 
CDT triangulation. 

These dimer models are of course somewhat artificial viewed from the per- 
spective of the two-dimensional CDT random triangulation. The rule for putting 
down the dimers (apart from being hard) is that you are not allowed to put them 
down on spatial links and for the links pointing forward in time from a vertex 
you are not allowed to put the dimer on the link furthest to the left. Nevertheless 
one would expect the critical behavior of a fullfledge hard dimer model on the 
random CDT surfaces to have the same critical behavior. This is what we will 
show below using the matrix model representation of a generalized CDT model 
introduced in [13]. 



The "plain" CDT model can be obtained as limit of an ordinary matrix model in 
the following way: assume the following matrix potential 



where the linear part of the potential is just for convenience of the scaling limit 
(see [T3] for a discussion), and where <fi is a N x N matrix. We then consider the 
following partition function: 



Expanding the potential in powers of A and performing the Gaussian integration 
results in an ensemble of random surfaces obtained by gluing together triangles 
(corresponding to the cubic term) and "tadpoles" corresponding to the linear 
term. The logarithm of Z represents the sum over the connected surfaces, and 
taking the large N limit will single out the surfaces of spherical topology. We will 
be interested in calculating the disk-amplitude, i.e. planar triangulations with a 
boundary. The corresponding object is 



3 The multicritical CDT matrix models 



3.1 Plain CDT 



(15) 




(16) 



W(x) = htr^—), 
JM x — <p 



(17) 
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Figure 4: The selfconsistent loop equation to be solved. The coupling constant 
g s associated with the double link monitors the tendency for the creation of baby 



universes. 



where the expectation value is wrt the partition function (16). In the large iV 



limit the disk amplitude W(x) satisfies the so-called loop equation 

g s W{xf = V'(x)W(x) -Q(x), Q(x)= Cl x + c , (18) 

where V denotes the derivative with respect to x. The parameter x is related to 
a boundary cosmological constant //& by x = e^ b . The graphic representation of 
the equation is shown in Fig. |4} 
The solution is 



W(x) 



V'(x)- y/V'(x) 2 -4g s Q(x) 
2g s 



(19) 



and the constants Co and c\ are determined by the requirement that the W(x) 
has a single cut on the real axis and that W(x) = 1/x + 0(l/x 2 ) for x — > oo. 
Thus W(x) has the following form 



W(x) 



V'(x) + X(x - c)y/{x -b){x- a) 
2g s ' 



(20) 



In this formula a, b and c are functions of the coupling constants A and g s . The 
continuum limit where W(x) can be associated with 2d Euclidean quantum grav- 
ity is, for a fixed g s , the limit where c(A) — > b(X). In the neighborhood of this 
point, A c , one has the expansion 



A = A c (l - e 2 )A, 6(A) = 6(A C ) - eVA, c(A) = c(A c ) + -eVA, (21) 



where 



x 



x c + eX, 



x r 



c(A c ) = 6(A C ). 



(22) 
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The interpretation is that A is the continuum cosmo logical constant, X the con- 
tinuum boundary cosmo logical constant. In this limit the term V'(x) does not 
scale and one obtains 

W(x) = NS(x) + const. e 3/2 W cont (X) , (23) 



W cont (X) = (X - y/K/2)\jX + y/X, (24) 

where W cont (X) agrees with the disk amplitude calculated using quantum Liou- 
ville theory and NS(x) is a non-scaling part, analytic in x. 

It is however possible to take another scaling limit related to CDT when 
g s — > 0, more specifically we have to scale g s as follows 

9s = G s e 3 . (25) 

In this limit, which was denoted the "classical" limit, the behavior of the critical 
couplings to lowest order in e is 

Uds) = \- \gT e\ x c {g s ) = 1 + G 1 / 3 e, (26) 

and the continuum cosmological constant and boundary cosmological constant 
are defined as 

A = \ c (g s ) - e 2 A, x = x c {g s )+eX. (27) 
Contrary to the situation for a fixed g s , W(x) itself will now scale and one obtains 

W(x)=e- 1 W cont (X), (28) 



WcoJX) = 2GT^ ' (29) 

where 

A cdt = A + ^G 2 J 3 , X cdt = X + G'J 3 , 2A cdt H - H 3 = 2G S . (30) 

We note for future reference that the critical points 

(A„ z*) := {X c {g s = 0),x c (g s = 0)) = (1/2, 1) (31) 
are characterized by 

V'(K,x m ) = V"(K,x*) = 0. (32) 
In the limit where G s — > we obtain precisely the CDT disk function: 
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X — > oo corresponds to contracting the boundary of the disk to a (marked) point, 
and the leading singularity of the corresponding closed (spherical) surface is thus 
y/X, or re-introducing the bare coupling constants, (A — A c (g s = 0)) 1 / 2 . This 
critical behavior is precisely the same as the critical behavior of the generic BP 
behavior, in agreement with the identification of the ensemble of CDT random 
surfaces with an ensemble of generic BPs mentioned above. 



3.2 Multicritical "classical" matrix models 

We now generalize the construction of a classical limit for "plain" CDT to the 
multicritical case. We consider the matrix potential 

V(<f>) = \<P 2 - A0 - ^ 3 - ^0 4 . (34) 

Viewed as a generating potential for random surfaces this matrix potential will 
glue together triangles and squares. Viewing each square as two triangles, we can 
think of the squares as part of the triangulation, but with a dimer placed on the 
diagonal, with a fugacity £ = A£. In this way the model describes dimers put on 
random triangulations in a special way, such that there is at most one dimer per 
triangle. On the graph dual to the triangulation the dimers are precisely hard 
dimers, and we will call them hard dimers also on the triangulation, even if the 
rule of putting down the dimers is slightly different from the standard hard dimer 
rule. Similarly we will denote £ the fugacity of the dimers, although it is strictly 
speaking £ which serves as the fugacity. 



We are interested in a multicritical behavior of the matrix model (34) and it 
occurs for negative fugacity £. The disk amplitude is still given by (19) where 
the polynomial Q(x) now is of second order with coefficients uniquely fixed by 
the requirement that we have a single cut on the real line and that W(x) falls off 
like \ jx for large x. The multicritical point X c (g s ), Cc(gs) is characterized by the 
following behavior of W(x), generalizing the critical behavior above: 

w( v V\x) - 2\l(g s )U9s)(x - 6 c (g.))V (s - c c (g s ))(x - a c (g s )) 

W{x) = — . (35) 

For a fixed g s the potential V (x) does not scale and plays no role in the continuum 
limit. 



We are interested in a "classical" limit of (35) where the potential scales, and 
that is obtained by the assignment 

g s = G s e\ (36) 



Note that the scaling in (36) differs from that in (25). This reflects that the 



multicritical point in the limit g s = is characterized by the generalization of 
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fl32l: 



V'{K,Z*,xJ = V"{K,^,x*) = V"'{K,^,x*) = 0. 



From (34) and (|35[) it follows that 



1 

A* 



2e* 



Vs. 



(37) 



(3? 



Thus a canonical scaling of the boundary cosmological constant around x* will 
lead to a V(A*, £*, x) of order £ 3 and a behavior compatible with (28 ) is obtained. 
With g s scaling like (36) one can calculate \ c (g s ),(, c (g s ),x c (g s ): 



A c (<? s ) = X (l - ^-G\l V - ^p^V) + 0(, 4 ) 

= 6 (i - fcl^ + ^G^e 3 ) + 0(. 4 ), 
x c (^) = W = (l + glrjG^e) + 0(e 2 ). 



(39) 



(40) 
(41) 



The perturbation away from X c (g s ), £ c (<?s) which leads to a potential V'(x, A, £) 
of order e 3 , assuming the boundary cosmological constant is perturbed as x = 
x c (g s ) + eX, can be parametrized as 



A = A* + As 2 - Ae 3 , £ = ^-^Ae 2 . 



(42) 



As in the ordinary multicritical matrix model situation one finds a two-parameter 
set of solutions depending on A, A, rather than the solution (29) which only 
depends on one paramter, A. Let us choose a convenient "background", using the 
notation from ordinary matrix models [IE], which we call the CDT-background, 
namely A = 0. By a redefinition similar to (30): 



A c dt — AH — G 3 / 4 , X c dt — X H — ;= — —GV 4 , 



>/3 5V4 



we obtain 



V\x) =£ 3 (A cdt + ^X, 



3 

cdt 



(43) 



(44) 



One can now calculate W(x) in the CDT limit G s — > where any creation of 
baby universes is suppressed: 



W(x) 



1 



1 



e *cdt + A& 8 ' 



(45) 
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Repeating the exercise below eq. (33), taking X cc ^ t to infinity, we obtain that the 
critical behavior of the partition function of spherical multicritical CDT surfaces 
with one marked point is dominated by the term (A — \ c (g s = 0)) 1//3 , in agreement 
with the expectation we have from the multicritical BP model. 

Finally, let us note that it follows from (42) that the critical exponent a of the 
hard dimers is 1/2, again in complete analogy with the multicritical BP model 
calculation (|6j)-( 12 ) . The result a = 1/2 is also the result obtained for ordinary 
multicritical matrix models. In this sense it is not really surprising that we obtain 
the result. One obtains the CDT multicritical point by following the ordinary 
multicritical line A c (<? s ), £, c {g s ) from, say, g s — 1 to the CDT value g s = 0. For any 
fixed value of g s we have of course o = 1/2, and although a might have changed 
when taking the particular classical scaling limit g s = e A G s , this is apparently 
not what happens. 



4 Discussion 

Only a few riddles are left in 2d Euclidean quantum gravity coupled to matter. 
One of them is the behavior of the Hausdorff dimension dh as a function of the 
central charge c of the conformal theory coupled gravity. A formula was derived 
by Watabiki some years ago [19] 



4 = 2^±^ (46) 

Most likely this formula is correct for c < 0. For c = agrees with what is 
known to be the correct answer [2D]. For c = —2 there are very reliable computer 
simulations which show agreement with the formula [21]. Finally for c — > — oo it 
gives 2, something one would indeed expect from semiclassical Liouville theory. 
However, for < c < 1 the numerical agreement is less conclusive [22], and the 
possibility that d h = 4 in this range was pointed out, and the idea has recently 
been resurrected [23] • For c > 1 the Liouville formulas become complex and 



expressions like (46) are not valid, but it is believed that there is a universal 
phase where the world sheet degenerates to BPs. 

Surprisingly we have a somewhat similar situation in CDT: from numerical 
simulations dh seems unchanged (and equal 2) when matter with central charge 
< c < 1 is coupled to the CDT ensemble [21] and recently it was shown that 
there might be a kind of universal phase for c > 1 |25j . However, to the extent 
we can really view the hard dimer models as corresponding to conformal field 
theories, it seems that for c < the matter systems can change fractal structure 
of the CDT ensemble. Qualitatively the changes are like in the full Euclidean 
models, d h decreases with decreasing c. In the c = case the CDT model can 
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be understood as an effective Euclidean model, where baby universes have been 
integrated out. Whether such an interpretation is possible also when matter is 
coupled to CDT is presently unknown, but since the multicritical matrix models 
capture the critical behavior of both CDT and ordinary 2d Euclidean quantum 
gravity coupled to certain matter systems, depending on how we scale g s , we have 
a chance to answer this question in the context of matrix models. 

Note 

While completing this article we were informed by Stefan Zohren that he and 
Max Atkin have obtained results which are identical to a number of our results. 
We thank Stefan for informing us of these results prior to publication. 
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